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Continuing  work  in  the  design  of  shock  ignition  targets  is  described.  Because  of  reduced  implosion  velocity  require¬ 
ments,  low  target  adiabats,  and  efficient  drive  by  short  wavelength  lasers,  these  targets  produce  high  gain  (>  100) 
at  laser  energies  well  below  1  megajoule.  Effects  of  hydrodynamic  instabilities  like  Rayleigh-Taylor  or  Richtmyer- 
Meshkov  are  greatly  reduced  in  these  low-aspect  ratio  targets.  Of  particular  interest  is  the  optimum  ratio  of  ignitor  to 
compression  pulse  energy.  A  simple  pellet  model  and  simulation-derived  coupling  coefficients  are  used  to  analyze  opti¬ 
mal  fuel  assembly,  and  determine  that  shock  ignition  allows  enough  control  to  create  theoretically  optimum  assemblies. 
The  effects  on  target  design  due  to  constraints  on  the  compression  and  ignitor  pulse  intensities  are  also  considered  and 
addressed.  Significant  sensitivity  is  observed  from  low-mode  perturbations  because  of  large  convergence  ratios,  but  a 
more  powerful  ignitor  can  mitigate  this. 


I.  INTRODUCTION 

Shock  ignition  is  a  new  concept  for  producing  inertial  con¬ 
finement  fusion  (ICF),  although  it  has  precursors  in  earlier 
research1.  Like  conventional  central  ignition  ICF,  it  relies  on 
a  small  central  hot  spot  that  ignites  and  burns  the  surrounding 
cold  high  density  fuel.  However,  in  shock  ignition  the  target 
compression  and  ignition  are  produced  by  two  distinct  parts  of 
the  laser  pulse:  a  compression  pulse  and  a  short,  high  intensity 
ignitor  pulse  respectively. 

Analysis  indicates  that  shock  ignition  can  produce  much 
higher  gains  at  low  laser  energies  than  conventional 
ignition2,3.  Simulations  find  that  significant  gains  could  be 
achieved  with  as  little  as  ~  250  kJ  of  0.248pm  wavelength 
light.  One -dimensional  (ID)  gains  of  100  are  possible  for 
laser  energies  of  about  a  quarter-megajoule,  and  gains  are 
comparable  to  fast  ignition  schemes4.  In  high  resolution  two- 
dimensional  (2D)  simulations  that  include  single  nonunifor¬ 
mity  sources,  targets  provide  nearly  full  ID  gains  with  ex¬ 
pected  levels  of  outer  perturbations,  inner  perturbations,  or 
ISI  imprint.  The  short  wavelength  and  zooming  capabilities  of 
the  Krypton  Fluoride  (RrF)  laser  result  in  much  more  efficient 
drive,  and  higher  gains,  than  frequency-tripled  glass  lasers. 

This  paper  summarizes  results  from  new  shock  ignition  de¬ 
sign  studies.  We  are  particularly  interested  here  in  the  optimal 
separation  of  the  laser  pulse  energy  into  the  target  compres¬ 
sion  and  hot-spot  ignition.  Like  fast  ignition,  shock  ignition 
promises  independent  control  of  heating  and  compression  of 
the  hot  spot  and  cold  fuel.  We  benchmark  a  simple  theoretical 
model  that  predicts  an  optimal  splitting  of  the  pellet  energy, 
and  investigate  the  coupling  between  the  compression  and  ig¬ 
nitor  pulses  and  the  pellet  assembly  components.  Such  control 
implies  that  there  is  a  continuum  of  designs  between  shock  ig¬ 
nition  and  conventional  central  ignition:  with  sufficiently  high 
compression  power  the  needed  ignitor  shock  vanishes  and  one 
recovers  conventional  central  ignition.  This  offers  much  flex- 
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ibility  in  the  design  of  high  performance  implosions. 

Because  of  the  reduced  implosion  velocity  requirements 
in  shock  ignition,  one  can  use  compression  intensities  that 
are  smaller  than  needed  for  conventional  ignition;  alterna¬ 
tively,  the  target  radius  and  aspect  ratio  (AR)  can  be  re¬ 
duced  to  produce  more  stability.  The  latter  approach  is  taken 
here,  so  intensities  are  limited  by  the  need  to  avoid  laser- 
plasma  instabilities  (LPI)  during  compression.  The  actual 
limits,  however,  are  uncertain.  Previous  studies  at  0.35 pm 
laser  wavelength  have  suggested  that  the  two  plasmon  de¬ 
cay  instability  begins  to  be  of  concern  in  the  neighborhood 
of  ~  5  —  8  x  1014  W/cm25~7 .  When  the  intensity  increases 
past  this  limit,  the  instabilities  show  nonlinear  behavior  and 
are  not  well  described  by  existing  theory.  Results  from  recent 
LPI  experiments  on  the  Nike  KrF  laser  (0.248 pm  wavelength) 
are  consistent  with  an  intensity  limit  of  about  a  factor  of  two 
higher8,9,  and  further  experiments  are  planned.  Whatever  the 
limits,  the  compression  intensity  determines  the  pressure  that 
drives  the  pellet,  and  ultimately  the  hydrodynamic  instability 
growth  and  yield  that  can  be  achieved. 

In  contrast,  the  high  shock  pressures  needed  for  ignition  re¬ 
quire  ignitor  intensities  higher  than  these  limits,  and  LPI  will 
likely  be  active.  When  the  ignitor  is  applied,  however,  the  tar¬ 
get  has  an  areal  density  large  enough  to  shield  itself  from  hot 
electrons  below  about  ~  100 keV10.  Indeed,  the  production  of 
hot  electrons  during  the  ignitor  may  actually  raise  the  ignitor 
coupling  efficiency11.  Thus  the  ignitor  may  be  constrained 
as  much  by  the  laser  capability  as  by  fundamental  interaction 
physics  issues.  In  this  paper,  we  do  not  simulate  LPI  effects 
but  instead  examine  the  tradeoffs  that  may  be  possible  to  avoid 
possible  intensity  constraints. 

The  results  presented  here  are  separated  into  two  parts;  the 
first  section  deals  with  design  work  in  ID  using  the  FAST  ra¬ 
diation  hydrocode12,  and  examines  limits  and  tradeoff  in  pa¬ 
rameters  determining  the  compression  and  ignitor  pulse  in¬ 
tensities.  Simulation  results  are  compared  to  a  simple  model 
of  the  assembled  pellet,  and  the  limits  of  constructing  opti¬ 
mal  configurations  are  examined.  Options  for  changing  the 
intensities  of  the  laser  pulses  are  also  examined.  The  second 
section  reports  on  2D  FAST  simulation  results,  and  considers 
the  effects  of  hydrodynamic  instabilities  and  limits  of  target 
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FIG.  1.  The  ID  gain  curve  constructed  from  the  optimum  gain  cases 
at  each  scale  factor  for  the  shock  ignition  targets  simulated  in  this 
study. 


perturbations.  Options  for  coping  with  the  high  convergence 
ratios  and  intensity  limitations  are  presented  and  discussed. 


II.  ID  STUDIES 

The  shock  ignition  targets  examined  here  are  DT-ice  fuel 
shells  within  low-density  wetted-and-frozen  DT  foam  abla¬ 
tors,  and  have  a  low  initial  aspect  ratio  ( R/AR  =  2.5).  The 
target  mass  is  varied  by  a  factor  of  16,  with  the  aspect  ratio 
and  initial  shock  pressure  remaining  the  same  (see  Table  I). 
A  relaxation  pulse13  drives  the  target.  An  initial  picket  con¬ 
taining  less  than  a  percent  of  the  total  laser  energy  induces  a 
~  50  Mbar  decaying  shock  through  the  target.  This  is  fol¬ 
lowed  by  a  foot  pulse  with  a  shock  pressure  of  ~  3  Mbar , 
driven  by  an  intensity  of  ~  6  TW /cm2.  Further  shock- free 
compression  leads  to  a  1/e  mass-averaged  adiabat  between 
1.8-2.  (The  adiabat  is  defined  as  the  ratio  of  the  pressure  to 
the  Fermi-degenerate  pressure).  This  shock-free  compression 
is  done  by  the  remaining  compression  pulse,  which  rises  to 
peak  powers  producing  pressures  of  100’s  of  Mbar.  Finally,  a 
short  pulse,  high  intensity  “spike”  is  applied  just  before  stag¬ 
nation,  producing  a  shock  of  order  1  Gbar.  The  ignitor  start 
is  contiguous  with  the  end  of  the  compression  pulse.  Laser- 
pellet  coupling  is  optimized  by  zooming  the  laser  (i.e.,  shrink¬ 
ing  the  focal  spot)  twice  during  the  implosion14. 

We  optimize  the  target  by  varying  the  compression  power 
and  the  ignitor4.  The  target  is  first  tuned  by  selecting  a  com¬ 
pression  power  (without  an  ignitor)  and  maximizing  the  peak 
areal  density.  Then,  a  variety  of  ignitor  pulses,  defined  by 
their  peak  power,  pulse  width,  and  turn-on  time  (“timing”) 
are  applied  and  gain  is  optimized.  The  compression  power  is 
changed  and  the  process  is  repeated.  This  maps  out  a  region 
in  a  four-parameter  space  that  produces  the  highest  gains.  Se¬ 
lecting  the  best  gain  achieved  at  each  pellet  scale  produces  the 
gain  curve  (Fig.  1). 


TABLE  1.  Target  specifications,  listed  by  relative  mass.  In  all  cases, 
the  ablator  is  made  of  100 mg /ctn?  CH  foam  wicked  with  solid  DT, 
the  fuel  layer  is  DT  ice.  The  target  linear  dimensions  are  scaled  by 
the  cube  root  of  the  scale  factor,  which  denotes  the  relative  target 
mass.  Simulation  parameters  listed  are  those  found  at  highest  gain. 


Scale  ( relative  mass ) 

lx 

2x 

5x 

8x 

16x 

Outer  Radius  (pm) 

854 

1076 

1455 

1708 

2152 

Ablator  Thickness  (pm) 

108 

136 

177 

216 

272 

Fuel  Thickness  (pm) 

237 

300 

407 

476 

600 

Best  gain 

101 

143 

210 

243 

309 

R laser  (kJ) 

231 

390 

727 

1074 

1700 

Absorption  (%) 

80 

85 

87 

91 

93 

IFAR 

21 

20 

16 

14 

11 

pR(g/cmi) 

2.4 

2.9 

3.4 

3.4 

4.0 

Ignitor  power  (TW) 

300 

400 

450 

500 

750 

Ignitor  pulse  width  (ps) 

200 

250 

500 

400 

550 

Compression  power  (TW) 

95 

110 

120 

160 

150 

In-flight  adiabat 

2.0 

2.0 

1.9 

1.7 

1.8 

k implosion  (*10  cm/s ) 

2.8 

2.6 

2.2 

2.2 

1.9 

R ignitor  /R\aser 

0.18 

0.19 

0.21 

0.13 

0.20 

Rr  \R()/ Rfiotspotmh,  \ 

63 

65 

68 

52 

65 

A.  Optimization  of  Compression  Power 

A  priori,  it  is  unknown  what  compression  power  produces 
the  best  gain.  The  implosion  velocity  should  be  smaller 
than  needed  for  ignition,  yet  high  enough  to  give  significant 
compression.  As  a  rule  of  thumb,  we  begin  using  compres¬ 
sion  powers  that  produce  implosion  velocities  of  ~  2  —  2.5  x 
107  cm  j sec  (for  the  scale-2  target,  the  minimum  ignition  ve¬ 
locity,  without  an  ignitor,  is  >  3  x  107  cm /sec). 

The  best  gain  typically  occurs  at  compression  powers  less 
than  50%  of  that  needed  for  conventional  ignition  (Fig.  2),  but 
the  optimized  gains  are  relatively  insensitive  to  changes  in  the 
compression  power.  As  the  power  increases,  the  compression 
pulse  alone  becomes  large  enough  to  ignite  the  pellet,  and  the 
ignitor  is  no  longer  needed.  The  gain  drops  towards  this  limit 
as  the  laser  energy  requirements  become  larger.  (This  ignitor¬ 
less  limit  would  be  precluded  in  practice  because  the  intensi¬ 
ties  needed  to  achieve  it  [Fig.  2b]  are  well  above  LPI  limita¬ 
tions.)  Thus,  there  is  a  continuum  between  shock  ignition  and 
conventional  central  ignition.  For  smaller  compression  pow¬ 
ers,  the  ignitor  energy  grows  and  the  timing  tolerance  narrows. 
As  the  ignitor  increases,  the  coronal  plasma  becomes  hotter, 
the  absorption  drops  and  the  shock  launch  time  moves  closer 
to  the  stagnation  time.  Simultaneously,  the  lower  compression 
power  produces  smaller  pR  and  lower  gains  follow. 

Predictions  of  optimum  compression  power  require  knowl¬ 
edge  of  (a)  the  optimum  configuration  of  the  imploded  and 
assembled  fuel  and  (b)  the  coupling  between  the  laser  pulse 
components  and  the  assembled  pellet  components.  We  con¬ 
sider  coupling  in  section  II B;  here  we  turn  our  attention  to 
the  optimum  fuel  assembly.  A  simple  Meyer-ter-Vehn-type 
analysis15  can  be  used  to  show  that  for  optimally  assembled 
targets,  the  hot  spot  energy  should  be  much  smaller  than  the 
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FIG.  2.  (Color  online)  (a)  Optimum  gain,  (b)  peak  laser  intensities, 
and  (c)  laser  pulse  energies  as  a  function  of  compression  power,  for 
a  250  ps  FWHM  ignitor  and  the  scale-2  target. 


compression  energy  (see  Appendix  A),  about  27  times  smaller 
for  an  isobaric  target.  This  hot-spot  to  cold-fuel  energy  ratio  is 
larger  for  shock  ignition  targets,  depending  upon  the  pressure 
ratio  between  the  hot  spot  and  the  cold  fuel. 

This  model  was  verified  with  simulations  of  different  scale 
targets  and  compression  powers,  using  drive  parameters  that 
are  optimal  for  the  corresponding  scale  size  and  compression 
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FIG.  3.  Comparison  of  simulation  and  predicted  hot-spot  (open  di¬ 
amonds)  and  cold-fuel  (solid  squares)  energies.  The  predicted  ener¬ 
gies  use  the  formula  from  Table  II  in  Appendix  A,  and  the  assembled 
fuel  masses,  stagnation  adiabats,  and  hot-to-cold  pressure  ratios  mea¬ 
sured  in  in  the  simulations. 


power.  To  separate  the  effects  of  the  pulse  heating  from  alpha 
particle  heating,  thermonuclear  burn  is  turned  off.  The  pellet 
is  considered  fully  assembled  when  the  internal  energy  of  the 
hot  spot  peaks.  (The  hot  spot  is  defined  as  the  fuel  around 
r  =  0  with  a  temperature  of  at  least  1  keV,  while  the  cold  fuel 
is  defined  as  the  mass  within  1  je  of  the  peak  density).  Fig.  3 
shows  that  the  hot  spot  and  cold  fuel  energies  from  the  simu¬ 
lation  scale  according  to  the  predictions  from  this  model. 

The  laser  compression  and  ignitor  pulses  separately  affect 
the  components  of  an  optimal  configuration.  In  the  next  sec¬ 
tion  we  examine  the  contributions  of  each  part  of  the  laser 
pulse  (compression  vs.  ignitor)  to  the  two  assembled  fuel  re¬ 
gions. 


B.  Coupling 

In  conventional  ignition,  the  laser  pulse  couples  a  fraction 
of  its  energy  to  the  pellet,  known  as  the  coupling  efficiency. 
However,  knowledge  of  this  single  fraction  is  insufficient  to 
optimize  the  pellet,  since  we  wish  to  independently  control 
the  pellet’s  hot  spot  and  cold  fuel  with  the  compression  and 
ignitor  pulses.  Both  of  these  two  pulses  produce  pellet  com¬ 
pression  and  heating,  but  in  different  ways.  The  compression 
pulse  compresses  and  accelerates  the  fuel  shell,  producing  ki¬ 
netic  energy  which  is  then  converted  into  internal  energy  when 
the  pellet  stagnates.  The  ignitor  delivers  energy  through  the 
spherically  convergent  shock  it  producing  shock  heating  and 
p  —  dV  work.  Because  of  their  very  different  intensities,  the 
absorption  is  also  markedly  different  for  the  two  pulses.  We 
wish  to  know  if  the  coupling  of  these  two  pulses  to  the  hot 
spot  and  fuel  is  different  enough  to  produce  a  truly  optimal 
pellet  configuration. 

We  first  construct  the  coupled  equations  that  relate  hot-spot 
and  cold-fuel  energy  to  the  compression  and  ignitor  pulse  en- 
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FIG.  4.  Absorption  and  coupling  efficiencies  for  the  laser  pulse  to  the  pellet  ignition  assembly.  The  coupling  efficiencies  have  the  absorption 
factored  out,  i.e.,  they  denote  coupling  of  absorbed  energies.  The  diamond  symbols  denote  the  laser  compression  pulse,  and  the  square  symbols 
denote  ignitor  pulse  coupling. 


ergies: 


&  pellet  —  Me  laser 


(i) 


where  C  pellet  —  hotspot  s^coldfuel\s  £ laser  —  \r  ignitors  ^compression] 

and 


M  = 


£/ii  &hc 
£ci  £cc 


(2) 


is  the  coupling  matrix.  If  the  matrix  is  diagonally  dominant, 
then  we  have  independent  control  of  the  energy  in  the  pellet 
hot  spot  and  the  cold  fuel. 

Next  the  laser  pulse  energies  are  found  as  eiaser  = 

=-i_ 

M  &  pellet  5  t)r  . 

^ ignitor  ( hot  spot  ^hc^coldfuel )  /  dctl.V/j  (3) 

C compression  (  ^ci^hotspot  T  ^hi^coldfuel  j  /  detl  /V/  j  (4) 

where  det(M)  =  £/„ece  —  £/)e£e;.  Thus  the  ratio  of  the  ignitor 
energy  to  the  laser  compression  energy  is: 


€  ignitor 
& compression 


£cc  eh:f  £hc 


(5) 


from  nascent  thermonuclear  burn,  burn  is  turned  off.  Simula¬ 
tions  with  and  without  ignitors  determine  its  additional  heat¬ 
ing.  Typical  values  (see  Fig.  4)  for  the  coupling  efficien¬ 
cies  show  that  the  coupling  matrix  is  diagonally  dominant: 
£/„■  ~  1  —  1 .5%  and  £/ic  ~  0.2  —  0.5%,  while  ecc  ~  5  —  6%  and 
£ci  ~  2  —  4.5%.  Although  the  coupling  is  not  completely  or¬ 
thogonal,  the  ignitor  preferentially  couples  its  energy  to  the 
hot  spot  while  the  compression  pulse  couples  more  to  the  cold 
fuel. 

The  limits  of  the  ratio  ehs/ecf  that  can  be  produced  by 
shock  ignition  can  be  approximated  by  using  the  following 
average  values  for  the  net  coupling  coefficients:  £/„■  =  0.01, 
£hc  =  0.003,  £cc  =  0.055  and  £cl  =  0.025.  This  predicts  that 
we  can  create  pellets  with  0.055  <  ehs/ecf  <  0.4.  Accord¬ 
ing  to  the  simple  model  in  appendix  A,  the  optimum  ratio 
of  ehs/ecf  should  be  «  <t>/27.  In  the  targets  simulated  here 
<J>  «  2  —  3  (without  burn),  so  we  expect  optimal  assemblies  to 
have  ehs/ecf  «  0.09  —  0.13.  If  the  coupling  coefficients  are 
assumed  to  be  roughly  constant,  we  would  conclude  that  the 
optimum  ignitor/compression  energy  ratio  would  be  between 
0.13  and  0.43.  The  observed  ratios  of  ~  0.25  (Table  I)  are  in 
this  range. 


C.  Laser  Intensity 


Assuming  all  the  coupling  coefficients  are  positive,  this  rela¬ 
tion  is  only  positive  in  the  range: 

£/zc/ fcVr  A  ?hs / ecf  ^  Ehi  / Eci  •  OV  (6) 

Shi/ £rz  ?hs  f  ?cf  £}i<: / £cc-  (2) 

For  a  diagonally  dominant  coupling  matrix,  only  the  first  con¬ 
dition  applies.  This  inequality  shows  that  there  is  a  limited 
range  of  hot-spot  to  cold-fuel  energies  in  the  pellet  assemblies 
that  can  be  formed. 

We  find  the  coupling  values  with  a  set  of  simulations  at 
each  pellet  scale.  Compression  powers  are  sampled  at  opti¬ 
mal,  low,  and  high  values.  To  separate  pulse  heating  effects 


There  is  an  upper  limit  on  drive  intensity  which,  if  ex¬ 
ceeded,  can  produce  destructive  laser-plasma  interactions 
(LPI).  This  limit  will  be  determined  by  the  plasma  conditions 
(primarily,  temperature  and  density  scalelength)  and  laser 
wavelength.  Although  the  exact  limit  is  currently  unknown, 
current  measurements  suggest  that  it  is  less  than  1015  W /cm2 
for  0.35  jlm  wavelength  light5-7,  and  a  factor  of  about  two 
higher  (1.5  —  1.7  x  1015  W /cm2)  for  0.248  jlm  wavelength 
light9.  The  primary  LPI  concern  in  direct-drive  is  the  two- 
plasmon  decay  instability  (TPD),  because  the  threshold  is  low 
and  it  can  generate  electrons  hot  enough  to  penetrate  and  pre¬ 
heat  the  fuel.  Stimulated  Raman  scattering  (SRS),  stimulated 
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Brillouin  scattering  (SBS)  and  the  parametric  decay  instability 
(PD)  are  possible  as  well. 

The  intensity  limit  is  most  important  during  compression, 
since  the  pellet  is  most  sensitive  to  preheat  during  this  time. 
The  compression  intensity  can  be  decreased  by  adjusting  the 
pellet  mass  and/or  aspect  ratio,  and  also  by  changing  the  tim¬ 
ing  and  zooming  ratios.  The  drawback  to  these  adjustments  is 
their  negative  impact  on  hydro  stability. 

However,  the  ignitor  pulse  intensity  must  climb  to  ~ 
1016  W  /cm2  to  produce  an  adequately  strong  igniting  shock. 
This  is  clearly  above  projected  LPI  thresholds.  As  has  been 
pointed  out  previously1011,  the  ignitor  begins  when  the  pel¬ 
let  areal  mass  has  significantly  increased  due  to  convergence. 
The  areal  mass  then  (~  0.05  —  0.06  gm/cm2)  can  shield  the 
pellet  interior  from  electron  temperatures  of  ~  100  keV 10.  LPI 
electrons  lower  than  this  limit  may  actually  help  to  transport 
the  laser  energy  deeper  into  the  pellet  ablator  and  enhance  the 
igniting  shock.  However,  if  other  LPI  effects  are  destructive 
(e.g.,  reflecting  too  much  light  from  absorption  regions,  or 
producing  extremely  hot  electrons),  it  may  be  necessary  to 
limit  the  laser  intensity. 

We  examine  the  flexibility  of  controlling  the  intensity  in 
these  different  parts  of  the  pulse  in  the  next  two  sections. 


FIG.  5.  (Color  online)  The  ignitor  shock  pressures  as  a  function  of 
time,  as  the  ignitor  pulse  width  changes.  The  ignitor  peak  power 
is  450  TW.  The  solid  lines  are  the  peak  shock  pressures,  while  the 
dotted  lines  are  the  pressure  at  the  laser  critical  surface.  Inset:  The 
pellet  yield  (black  diamond)  and  gain  (red  triangle)  vary  with  the 
ignitor  pulse  width. 


1.  Ignitor  pulse  flexibility 

The  ignitor  pulse  uses  high  intensities  in  order  to  produce 
the  igniting  shock.  It  can  produce  LPI  but  also  heats  the 
plasma  corona,  which  decreases  collisional  absorption.  There 
may  also  be  driver  limitations  on  either  the  ignitor  pulse  length 
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FIG.  6.  The  gain  vs.  ignitor  timing  and  power  (and  pulse  length), 
holding  the  ignitor  energy  fixed  at  90  kJ. 


or  power16.  Longer,  lower  power  ignitors  may  be  easier  to 
generate,  but  short  higher  power  ignitors  may  couple  better. 
Here  we  examine  the  flexibility  offered  by  trading  off  pulse 
power  for  pulse  length  at  constant  ignitor  energy. 

After  the  ignitor  pulse  launches  its  shock,  it  turns  off  and 
a  rarefaction  wave  follows  the  shock.  Geometric  convergence 
increases  the  shock  strength  and  accelerates  it  as  it  nears  the 
target  center.  The  rarefaction  wave  chases  the  shock  and  can 
overtake  it  and  reduce  its  strength.  For  short  pulse  lengths  the 
ignitor  shock  will  not  have  time  to  reach  full  strength,  and/or 
the  rarefaction  wave  will  catch  and  erode  it.  Increasing  pulse 
lengths,  after  attaining  ignition,  can  increase  the  pellet  areal 
mass  and  increase  burnup,  but  expend  more  energy  and  can 
reduce  the  gain.  Thus,  there  may  be  some  optimum  pulse 
width. 

First,  we  varied  the  ignitor  pulse  width  at  constant  power 
(Fig.  5)  in  an  optimally  timed  scale-2  target.  The  ignitor 
pulse  rises  linearly  for  100  psec,  and  is  then  shut  off  at  vary¬ 
ing  times.  While  the  pressure  at  the  critical  surface  responds 
immediately  to  the  laser  pulse  turn-off,  the  peak  shock  pres¬ 
sure  is  unaffected  once  the  pulse  width  exceeds  120  psec.  By 
delaying  the  effects  of  the  incoming  rarefaction  wave,  larger 
pulse  widths  produce  more  fuel  compression  and  areal  mass, 
which  increases  the  yield.  However,  the  higher  yield  is  more 
than  offset  by  the  increase  in  laser  energy  used  in  extending 
the  ignitor.  At  constant  power,  the  optimum  width  is  just  a 
little  larger  than  the  minimum  needed  to  induce  ignition.  Al¬ 
though  it  is  not  evident  here,  a  larger  pulse  width  can  be  help¬ 
ful  in  increasing  the  pellet  robustness  to  2D  perturbations. 

There  is  also  a  broad  flexibility  in  trading  off  the  ignitor 
power  with  its  pulse  width,  at  constant  energy.  Using  the  same 
target,  we  fix  the  energy  at  90  kJ  (approximately  20%  of  the 
compression  energy,  and  typically  the  amount  producing  the 
optimum  gain).  The  results  (Fig.  6)  show  that  the  pulse  length 
can  be  varied  from  450  ns  to  less  than  ~  4  ps  while  still  keep¬ 
ing  effectively  the  same  gain;  for  most  of  these  powers,  the 
ignitor  timing  has  a  350  psec  range.  At  the  high  power/short 
pulse  width  end  of  this  parameter  range,  the  physics  as  simu- 
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lated  here  appears  to  be  sensitive  only  to  the  deposited  energy. 
In  reality,  such  high  powers  would  increasingly  be  affected  by 
LPI,  which  is  beyond  the  scope  of  this  modelling. 


2.  Compression  intensity 

Because  the  incoming  velocity  requirements  for  shock  igni¬ 
tion  are  significantly  lower  than  those  in  conventional  ignition, 
there  is  possibility  of  using  lower  compression  intensity.  On 
the  other  hand,  for  a  given  implosion  velocity,  a  higher  com¬ 
pression  intensity  will  allow  one  to  use  more  mass  in  the  tar¬ 
get  (producing  higher  pellet  yields)  and/or  to  use  lower  aspect 
ratio  targets  (producing  more  hydrodynamic  stability).  Thus, 
there  is  still  a  strong  incentive  to  use  the  highest  intensity  pos¬ 
sible.  In  the  targets  presented  so  far,  we  have  used  compres¬ 
sion  intensities  of  about  2  x  1015W / cm2 .  Here  we  investigate 
the  tradeoffs  involved  in  lowering  this  intensity  with  the  scale- 
2  target. 

The  target  was  rescaled  using  the  relations  R  ~  /  1/7  and 
pAR  ~  /°-45, 17  which  assume  constant  energy  and  velocity. 
After  tuning  the  new  targets,  the  ignitor  is  optimized  for  gain. 


Aspect  Ratio 

FIG.  7.  (Color  online)  The  maximum  laser  intensities  and  gains  for 
scale-2  pellets  rescaled  to  different  initial  aspect  ratios.  All  targets 
have  the  same  laser  energy  and  implosion  velocity  (without  ignitor). 

The  lowered  intensity  leads  to  a  larger  radius  target  with 
a  larger  initial  aspect  ratio:  to  achieve  the  same  implosion 
velocity  (about  2.4  x  107  cm/ sec),  the  pellet  shell  must  be 
pushed  further  with  the  decreased  pressure.  Although  the  shell 
is  thinner,  it  will  have  more  mass  (m  ~  R2AR  ~  /  >l22).  Lower 
intensity  also  results  in  cooler  temperatures  and  smaller  ab¬ 
lation  velocities,  which  give  higher  absorption  and  hydrody¬ 
namic  efficiency.  All  of  this  leads  to  higher  ID  yield  and  gain 
(Fig.  7).  The  drawback  is  that  the  pellet  shell  becomes  more 
susceptible  to  hydrodynamic  instabilities:  it  has  more  insta¬ 
bility  growth  as  the  shell  is  driven  longer  distances,  and  the 
smaller  mass  ablation  rate  leads  to  less  stabilization  for  high 
frequency  perturbations.  We  shall  return  to  this  point  in  the 


2D  analysis. 

III.  2D  STUDIES 

Two  dimensional  simulations  are  needed  to  assess  the  ef¬ 
fect  of  nonuniformity  sources  and  stability  of  the  targets. 
The  nonuniformities  include  pellet  manufacturing  imperfec¬ 
tions  (“surface  perturbations”)  as  well  as  driver  sources  such 
as  laser  imprint  or  beam  misalignment.  These  perturbations 
serve  as  seeds  for  hydrodynamic  instabilities  like  Rayleigh- 
Taylor  and  Richtmyer-Meshkov. 

Both  low  and  high  resolution  two  dimensional  simulations 
are  useful  here.  Low  resolution  simulations  (typically  850  ra¬ 
dial  zones  and  64  zones  in  inclination  between  0  =  0-  n) 
are  used  to  simulate  low-mode  nonuniformities;  the  range 
l  =  1  —  16  is  then  resolved  (8  cells  are  needed  to  resolve 
a  wavelength).  As  these  are  less  computer  intensive,  more 
simulations  are  available  to  determine  parameters  or  explore 
ranges.  High-resolution  simulations  (850  radial  by  2048  theta 
zones,  resolving  £=  1  —  5 12)  are  used  to  examine  small  wave¬ 
length  modes.  They  are  needed  to  resolve  growth  of  hydro- 
dynamic  instabilities,  which  can  grow  faster  at  higher  spatial 
frequencies. 


FIG.  8.  (Color  online)  Relative  rms  perturbation  of  the  areal  mass 
( SpR/pR )  in  the  pellet  for  the  scale-2  target  with  a  300  kj,  110  TW 
compression  pulse.  The  solid  black  line  shows  the  growth  of  outer 
surface  perturbations  (initially  0.48pm  rms),  the  red  dash-dotted  line 
is  the  inner  surface  perturbation  (initially  1  pm  rms),  the  blue  dash 
line  is  for  300  overlapped  1  THz  bandwidth  ISI  beams,  all  for  a  150 
kj  ignitor.  The  green  line  is  for  all  three  perturbation  sources  simul¬ 
taneously,  and  with  a  200  kJ  ignitor.  The  vertical  gray  lines  denote 
the  times  (a)  as  the  compression  shock  breaks  out  of  the  rear  of  the 
fuel  shell;  (b)  as  the  ignitor  shock  turns  on,  and  (c)  when  the  pellet 
stagnates  and  ignition  occurs. 

Nominal  surface  perturbations  (both  inner  and  outer)  and 
laser  optical  smoothing  (ISI)  generated  perturbations  are  on 
scale-2  targets  in  these  studies.  These  nominal  levels  are  as 
smooth  as  achievable  in  the  real  world.  The  outer  surface  per¬ 
turbations  correspond  to  the  surfaces  that  has  been  achieved 
and  measured  in  GDP  (glow  discharge  polymer)  plastic 
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FIG.  9.  (Color  online)  The  areal  mass  perturbation  amplitude  spectrum,  normalized  to  the  average  areal  mass  (A(f  =  0))  at  four  times  during 
the  implosion,  for  the  different  perturbations  sources  included  here.  The  different  times  shown  are  at  (a)  t=0;  (b)  12.8  ns  when  the  shocks 
break  out  from  the  rear  of  the  target;  (c)  16  ns,  when  the  ignitor  shock  turns  on;  and  (d)  16.9  ns,  when  the  pellet  ignites.  The  sources  included 
are  perturbation  on  the  outer  surface  (black),  inner  surface  (red  dash-dotted),  ISI  laser  imprint  (blue  dashed),  and  all  three  sources  combined 
at  once  (green).  The  modes  greater  than  t  =  10  have  been  locally  averaged  in  wavenumber  from  0.8f  to  1.2f  to  reduce  statistical  fluctuations 
that  obscure  the  plot. 


shells18,  and  have  an  amplitude  given  by  the  formula  A(\nm\  = 
9600£-3'39  +  2.4£~0-78  for  t  >  5  and  A([nm]  =  min(  102  — 
4.8(T  -  1.8)2'4, 100)  for  i  <  5.  This  is  a  0.125  jlm  {I  = 
2  —  2000)  rms  perturbation  in  plastic  (~  1  —  1.2  g/cm3), 
or  an  equivalent  areal  mass  perturbation  of  0.48  jlm  in  the 
0.33  gm/cm 3  wicked-foam.  The  inner  DT-ice  surface  per¬ 
turbation  has  the  form  A(\nm\  =  1000. /(3. Of0  6  +  2. 2-7f*)19, 
and  a  nominal  amplitude  of  1  jlm  rms.  For  optical  smooth¬ 
ing,  we  use  1  THz  ISI,  appropriate  for  KrF  light.  (This  is 
similar  to  1  THz  2D-SSD  smoothing  in  glass  lasers,  with 
the  exception  that  the  longer  wavelength  modes  (t  <  20)  are 
much  more  smoothed  with  ISI20).  300  overlapped  incoherent 
beams  are  chosen  to  correspond  to  NRL’s  Fusion  Test  Facility 
proposal21,  which  has  1800  total  beams.  (The  rule  of  thumb 
being  used  is  that  the  average  intensity  at  any  given  spot  on  an 
evenly  illuminated  pellet  is  N/6  of  the  intensity  from  a  single 
beam22,  where  N  is  the  total  number  of  beams  illuminating 
the  pellet). 

The  simulations  predict  that  target  is  robust  to  the  separate 
perturbation  sources  at  nominal  amplitude  levels,  with  ignitor 
power  and  energy  levels  comparable  to  but  greater  than  used 
to  produce  the  highest  ID  gain.  For  instance,  the  ID  best  gain 
was  found  to  be  142  for  450  TW,  200  psec  full  width  at  half 
maximum  (FWHM)  ignitor.  2D  high-resolution  simulations 
with  nominal  outer  surface  perturbations  on  this  target  find  an 
optimum  gain  of  125  for  a  650  TW  200  psec  ignitor  (the  ID 


best  gain  at  650  TW  is  138). 

For  a  fixed  laser  pulse  with  a  150  kJ  ignitor,  the  inner  DT- 
ice  perturbations  can  be  as  large  as  3  jlm  rms  before  ignition 
fails.  An  outer  surface  perturbation  at  twice  the  nominal  level 
will  spoil  ignition.  However,  doubling  the  ignitor  energy  to 
300  kJ  recovers  ignition,  as  the  gain  drops  from  120  to  94 
mostly  because  of  this  added  energy.  In  general,  as  the  pertur¬ 
bation  levels  increase  ignition  eventually  fails  as  the  hot  spot 
becomes  too  perturbed  and  insufficiently  hot  to  ignite.  For 
perturbation  levels  that  are  not  too  large,  increasing  the  igni¬ 
tor  power  restores  ignition  and  gain. 

Because  of  the  low  initial  and  in-flight  aspect  ratio  (IFAR), 
the  target  is  relatively  insensitive  to  the  usual  high-mode 
Rayleigh-Taylor  growth.  Separate  simulations  including  nom¬ 
inal  outer  or  inner  perturbations  or  ISI  imprint  show  that  mode 
growth  is  not  large  enough  to  penetrate  the  pellet  shell,  even 
when  amplified  by  the  ignitor  shock  (Fig.  8).  The  rms  per¬ 
turbation  level  of  these  individual  components  is  below  10% 
until  stagnation  and  burn  occur.  All  experience  a  brief  jump 
in  their  amplitude  when  the  ignitor  shock  is  applied,  and  gen¬ 
erally  more  rapid  growth  when  the  shell  decelerates  at  stagna¬ 
tion.  However,  the  laser  imprint  evolution  is  notable  because 
it  peaks  immediately  after  the  ignitor  shock  is  applied,  and 
then  monotonically  declines  as  the  shell  stagnates. 

Differences  in  the  perturbation  growth  spectra  can  explain 
this  divergent  behavior  (Fig.  9).  At  t=0,  there  is  no  per- 
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FIG.  10.  Images  of  an  imploding  pellet  at  times  (a)  100  psec  after  the  ignitor  starts;  (b)  as  the  ignitor  shock  is  halfway  through  the  shell; 
(c)  at  stagnation  as  the  burn  begins.  The  pellet  is  initially  perturbed  on  the  outer  and  inner  surfaces  with  “NIF-spec”  spectra4  with  nominal 
amplitudes  of  0.48  jlm  and  1.0  jlm  respectively,  and  was  subject  to  laser  imprinting  from  300  overlapped  mutually  incoherent  1  THz  ISI 
beams.  The  resulting  gain  is  102. 


turbation  yet  in  the  laser  imprint  case;  meanwhile  the  in¬ 
ner  and  outer  perturbation  levels  are  roughly  equivalent  for 
£  <  10  while  the  inner  perturbation  level  is  much  larger  at 
higher  modes.  At  shock  breakout  (near  the  end  of  Richtmyer- 
Meshkov  growth  and  before  Rayleigh-Taylor  (RT)  instabil¬ 
ity  begins),  both  the  outer  surface  and  imprint  perturbations 
have  grown,  while  the  inner  perturbation  is  unchanged.  The 
imprint  has  grown  larger  than  the  other  sources  in  the  range 
£  ~  10  —  200. 

When  the  ignitor  turns  on,  the  pellet  has  been  driven  in¬ 
ward  to  approximately  40%  of  its  original  radius.  The  growth 
in  this  interval  is  driven  mostly  by  the  ablative  RT  instability 
in  combination  with  Bell-Plesset  effects23  (Fig.  9c).  Imprint 
is  still  dominant  in  the  range  £  ~  10  —  100,  while  the  outer  per¬ 
turbation  dominates  for  smaller  and  larger  modes.  The  inner 
surface  perturbation  has  grows  for  £  ~  5  —  70,  but  decays  for 
lower  modes.  This  is  because  growth  of  the  amplitude  there 
is  smaller  than  the  growth  in  pR  due  to  the  spherical  conver¬ 
gence. 

Finally,  at  stagnation  (Fig.  9d),  all  growth  that  affects  the 
ignition  and  yield  has  occurred.  Compared  to  the  ignitor  on¬ 
set  some  900  psec  earlier,  the  spectra  are  marked  by  growth  in 
low  modes  (£  <  10  —  30)  and  decay  in  higher  modes.  This  is 
consistent  with  the  picture  of  deceleration  RT  growth  on  the 
inner  pellet  surface24:  higher  modes  are  strongly  stabilized  by 
the  ablation  driven  by  the  forming  hot  spot.  It  is  also  affected 
by  the  smaller  feedthrough  that  the  higher  modes  experience 
between  the  outer  shell  (where  they  were  driven  by  the  accel¬ 
eration  stage  RT)  to  the  inner  edge. 

Simulations  with  the  nominal  outer,  inner,  and  laser  imprint 
sources  combined  find  that  the  target  will  ignite  and  burn  well, 
but  only  if  the  ignitor  energy  is  increased  further.  For  example, 
the  simulation  shown  in  Fig.  10  gives  gain  ~  100  for  a  1  PW, 
200  psec  ignitor.  Again,  the  dominant  mode  at  the  end  of  the 
implosion  is  the  long  wavelength  £  =  5.  The  modal  spectrum 
of  this  simulation  is  also  shown  in  the  Figs.  9.  This  spectrum 
roughly  forms  an  envelope  over  all  the  other  individual  per¬ 
turbations,  and  confirms  that  the  outer  surface  perturbation  is 


the  main  source  of  the  final  low  mode  structure. 


relative  perturbation  amplitude 


FIG.  1 1 .  (Color  online)  Low-mode  (£  =  2  —  8)  simulations  for  target 
with  two  different  convergence  ratios.  The  perturbation  amplitude  is 
measured  as  a  ratio  to  the  nominal  perturbation  level  (0.48  pm  rms  in 
the  wetted  foam).  Solid  lines  are  the  laser  energy  used  at  the  optimal 
gain  (left  axis),  while  the  dotted  lines  give  the  optimal  gain  (right 
axis). 


A.  Convergence  Ratio 

While  the  targets  here  are  more  stable  due  to  their  low 
IFAR,  their  convergence  ratio  (Cr)25  is  relatively  high  (Ta¬ 
ble  I).  High  C,-targets  are  more  sensitive  to  low  mode  tar¬ 
get  perturbations,  which  is  consistent  with  the  dominant  low 
mode  geometry  in  the  stagnated  density  images  (Fig.  10c). 
These  long  wavelength  modes  grow  relatively  little:  classical 
RT  growth  is  limited  to  <  \J t  e-foldings,  about  a  factor  of  9 
for  £  =  5.  Other  growth  due  to  Richtmyer-Meshkov26  or  Bell- 
Plesset23  is  also  expected  to  be  small.  However,  even  if  the 
perturbation  remains  the  same,  the  relative  level  increases  as 
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~  1  /  R.  Pressure  and  mass  differences  lead  to  additional  sec¬ 
ular  growth,  and  scale  with  C,  (see  appendix  B).  These  con¬ 
siderations  predict  that  the  initial  relative  areal  mass  nonuni¬ 
formity  of  a  directly-driven  pellet  should  be  less  than  about 
0.1  /Cr;  for  the  observed  range  of  Cr  =  50  —  65  this  trans¬ 
lates  to  dpR/pR  <  0.15  —  0.20%.  The  outer  surface  finish 
of  0.48  pm  rms  in  the  ablator  here  corresponds  to  a  0.13% 
relative  areal  mass  perturbation,  and  is  close  to  this  predicted 
limit. 

We  try  to  mitigate  this  sensitivity  by  redesigning  the  scale-2 
target  to  have  a  smaller  Cr.  The  approach  taken  here  was  to 
re-time  the  shocks  from  the  foot  pulse  and  the  main  pulse  to 
meet  slightly  inside  the  fuel  shell  inner  edge.  This  increases 
the  adiabat  there,  and  pushes  more  fuel  into  the  low  density 
interior  that  later  forms  the  hot  spot.  Advancing  the  timing 
of  the  peak  compression  power  by  about  1 .2  nsec  decreased 
C,  from  65  to  about  49  for  the  scale-2  target  with  110  TW 
compression  power.  It  simultaneously  decreased  the  highest 
ID  gain  observed,  from  ~  140  to  ~  100. 

Lower  resolution  (£  ~  2  —  8)  2D  simulations  investigate  the 
effect  of  this  redesign  on  the  stability,  and  the  optimum  gain  is 
found  as  a  function  of  perturbation  amplitude.  Fig.  1 1  shows 
the  results,  which  do  not  reveal  any  advantage  to  lowering  Cr 
by  this  method:  the  gain  is  always  lower,  even  at  larger  pertur¬ 
bations  values.  The  increase  in  laser  energy  needed  for  larger 
perturbations  at  both  convergence  ratios  highlights  the  pre¬ 
ferred  method  of  handling  perturbations:  increase  the  energy 
in  the  ignitor. 


FIG.  12.  (Color  online)  Growth  of  outer  surface  modes  from  high- 
resolution  simulations  comparing  the  original  lower  aspect  ratio  tar¬ 
get  (AR=2.5)  with  a  target  that  is  designed  for  about  40%  lower  drive 
pressure  (AR=3.74).  The  cross  marks  on  the  curves  denote  the  time 
of  ignitor  turn-on  and  ignition. 


B.  Aspect  Ratio 

In  section  II C  2,  we  investigated  different  targets  based  on 
lower  drive  pressure  and  intensity.  There,  we  observed  that 
lower  intensities  lead  to  higher  gains  in  ID,  but  the  resulting 


FIG.  13.  Images  of  the  imploded  density  just  as  the  pellet  begins  to 
burn  (gain  is  ~  0.5  at  these  times)  for  (a)  lower  aspect  ratio  target 
(AR=2.5)  and  (b)  a  higher  aspect  ratio  target  (AR=3.7). 


higher  aspect  ratio  targets  were  also  projected  to  be  less  hydro- 
dynamically  stable.  Here  we  quantify  this  greater  instability. 
High  resolution  (Jt  =  1  —  512)  is  used  resolve  the  RT  growth 
during  the  implosion  of  ~  500 kJ  target  with  a  nominally  per¬ 
turbed  outer  surface.  This  target  has  an  initial  aspect  ratio  of 
3.74,  and  is  driven  with  a  maximum  compression  intensity  of 
about  6.6  x  1014  W /cm2  of  0.248  pm  light.  (For  0.35  pm 
laser  light,  it  requires  7  x  1014  W /cm2). 

Comparing  this  to  the  original  lower  AR  target  (Fig.  12), 
the  relative  rms  perturbation  growth  here  is  roughly  twice  as 
large  from  the  start  of  the  ignitor  through  to  stagnation.  Also, 
the  mode  structure  of  the  perturbations  exhibits  an  obvious 
shift  to  smaller  wavelengths  (see  Fig.  13),  indicating  that  the 
Rayleigh-Taylor  instability  is  becoming  more  important.  To 
attain  optimum  gain,  the  ignitor  energy  must  be  substantially 
increased  to  overcome  this.  For  the  original  AR=2.5  target, 
the  minimum  ignitor  energy  was  increased  from  70  kJ  (ID) 
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to  130  kJ  (2D)  to  ignite  and  produce  at  least  gain  100.  For 
the  AR=3.74  pellet,  the  minimum  ignitor  energy  needed  to 
achieve  gain  >  100  had  to  be  increased  from  57  kJ  in  ID  to 
225  kJ  in  2D.  In  addition,  the  maximum  gain  observed  in  2D 
for  this  latter  case  was  about  110,  whereas  gains  of  125  were 
seen  for  the  lower  aspect  ratio  target  with  the  same  perturba¬ 
tions. 


IV.  CONCLUSIONS 

We  have  investigated  low-adiabat  shock  ignited  targets  de¬ 
signed  to  be  driven  by  short  wavelength  lasers.  We  use  targets 
imploded  with  intensities  that  are  close  to  the  limits  imposed 
by  LPI,  in  order  to  minimize  the  target  aspect  ratio  and  max¬ 
imize  the  hydrodynamic  stability.  We  find  a  robust  parame¬ 
ter  space  for  driving  these  targets  as  functions  of  compression 
power,  ignitor  power,  and  timing.  The  robustness  monotoni- 
cally  increases  as  the  energy  in  the  ignitor  is  increased,  both 
in  the  ID  and  2D  simulations.  The  efficiency  of  energy  trans¬ 
fer  from  the  component  parts  of  the  pulse  (the  compressor  and 
the  ignitor)  is  also  calculated  from  the  simulations.  Although 
the  coupling  is  relatively  low  for  both  pulse  components,  the 
ignitor  couples  much  more  of  its  energy  into  the  hot  spot  of 
the  assembled  pellet  than  the  compression  pulse.  Thus,  the 
rudimentary  control  needed  to  separately  heat  the  igniting  hot¬ 
spot  independently  of  producing  compression  in  the  cold  fuel 
is  present.  This  control  allows  one  to  produce  near  optimal 
pellet  assemblies,  and  to  ignite  pellets  at  lower  drive  energies 
than  possible  with  conventional  central  ignition.  As  a  result, 
we  find  much  higher  gains  for  laser  energies  of  less  than  a 
megajoule.  For  optimal  gain,  the  ignitor  pulse  takes  up  about 
20%  of  the  total  laser  pulse  energy. 

The  low  aspect  ratios  of  these  pellets  make  them  more  ro¬ 
bust  to  hydrodynamic  instabilities.  This  is  possible  because 
the  implosion  velocities  are  low  (of  order  2  x  107  cm/ sec ),  and 
the  pressure  applied  by  the  compressing  laser  is  high  (~  100  — 
300  Mbar ).  The  latter  can  only  be  attained  in  practice  if  the 
compression  laser  intensity  used  (~  1.5-2  x  1015  W /cm2) 
is  not  significantly  affected  by  LPI.  While  this  is  more  likely 
to  be  the  case  for  the  0.248  jlm  laser  wavelength  used  here, 
determining  the  limiting  laser  intensity  is  still  an  area  of  ac¬ 
tive  research.  Current  experimental  results  indicate  that  this 
wavelength  may  allow  one  to  use  an  intensity  that  is  a  factor 
of  ~  2  higher  than  the  commonly  used  0.35  jlm  frequency- 
tripled  glass  laser  wavelength.  This  work  will  hopefully  spur 
research  needed  to  determine  these  intensity  limits  more  ac¬ 
curately. 

In  these  targets,  hydrodynamic  instabilities  produce  only 
mild  degradations  of  the  implosion.  Instead,  low-mode  sym¬ 
metry  constraints  due  to  the  target’s  large  convergence  ratios 
appear  to  be  the  dominant  multidimensional  effect.  In  general, 
these  low-mode  asymmetries  narrow  the  successful  operating 
space  allowed  for  these  targets.  Decreasing  the  convergence 
ratio  through  shock  timing  is  not  found  to  be  a  successful 
strategy.  Instead,  a  more  powerful  ignitor  can  overcome  these 
flaws,  although  the  gain  declines  somewhat  due  to  the  added 
laser  energy. 


If  the  drive  pressure  must  be  reduced  because  of  LPI  im¬ 
posed  limitations  on  the  compression  intensity,  then  the  pellet 
must  have  a  higher  initial  aspect  ratio.  This  leads  to  more 
susceptibility  to  hydrodynamic  instability.  While  we  have  not 
done  an  in-depth  study  of  these  tradeoffs,  we  have  explored 
the  effect  of  reducing  the  intensity  by  about  40%.  We  find  that 
the  initial  aspect  ratio  then  needs  to  be  increased  by  ~  50%, 
and  the  outer  radius  by  ~  20%.  This  higher  aspect  ratio  tar¬ 
get  allows  more  mass  (for  fixed  areal  density  and  laser  en¬ 
ergy)  and  thus  more  gain  in  ID.  However,  the  target  is  more 
Rayleigh-Taylor  unstable  in  2D,  and  the  higher  ignitor  energy 
needed  to  ignite  the  more  perturbed  hot  spot  leads  to  a  reduced 
gain  overall. 
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Appendix  A:  Pellet  Assembly  Analysis 

We  follow  here  the  analysis  of  Meyer-ter-Vehn  and 
others15,27-28  and  approximate  the  assembled  fuel  as  a  homo¬ 
geneous  hot-spot  surrounded  by  a  homogenous  mass  of  com¬ 
pressed  fuel.  The  main  difference  in  the  analysis  presented 
here  and  earlier  ones  is  that  we  vary  the  pressure  of  the  com¬ 
pressed  fuel  instead  of  the  hot  spot  radius15  27  or  density28. 
We  explicitly  differentiate  between  the  minimum  energy  re¬ 
quired  for  ignition  and  the  maximum  gain,  and  concentrate 
on  the  partitioning  of  the  energy  between  the  hot  spot  and 
the  cold  fuel.  Finally,  we  will  examine  the  differences  be¬ 
tween  the  isochoric  (fast  ignition),  isobaric  (conventional  di¬ 
rect  drive)  and  “mixed”  (shock  ignition)  schemes. 

The  hot  spot  energy  and  mass  are  determined  by  the  two 
conditions  required  for  ignition:  (1)  the  ion  temperature  must 
be  greater  than  a  minimum  ignition  temperature,  7}on  >  Tjgn  ~ 
5 keV,  and  (2)  the  hot  spot  areal  density  must  be  large  enough 
to  trap  3.5 Mev  alpha  particles:  pR/,  >  0.3g/cm2.  Assuming 
that  the  pressure  and  energies  in  the  hot  spot  can  be  described 
by  the  ideal  gas  law,  and  the  ion  and  electron  temperatures  are 
equilibrated,  the  hot  spot  energy  is: 

Eh  [ kJ ]  =  9.54 x  10 5pR3T3ph2bar.  (Al) 

T  and  pR  are  the  quantities  7}0„  and  pR/,  normalized  to  the 
nominal  values  5  keV  and  0.3 g/cnr. 

The  internal  energy  in  the  cold  compressed  fuel  is  also  de¬ 
termined  by  its  pressure.  The  two  constraints  that  we  need, 
corresponding  to  the  Tign  and  pRh  in  the  hot  spot,  are  the 
cold  fuel  mass  Mp  and  the  fuel  adiabat  a  =  p/pfd  where 
ppo[Mbar\  =  2.2 p5/li  is  the  Fermi-degenerate  pressure  and 
£pD[MJ/g\  =  0.33p2  ’’  is  the  Fermi-degenerate  energy.  Then 

Ec  [ kJ }  =  3.77 a°-6Mfmgp°G4bar  (A2) 
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where  Mfmg  is  the  fuel  mass  in  mg.  a  includes  ion  pressure, 
and  is  the  stagnation  adiabat,  not  the  in-flight  adiabat29'30. 
Any  remaining  kinetic  energy  in  the  imploding  shell  is  ig¬ 
nored,  so  in  this  sense  this  is  an  analysis  of  the  “optimum” 
configuration.  Here  and  in  the  following,  the  fuel  mass  is  as¬ 
sumed  to  be  in  units  of  mg,  the  pressure  in  Gbar,  and  the 
energy  in  kJ. 

The  total  energy  of  the  assembled  pellet  is  the  sum  of  the 
hot-spot  and  cold-fuel  energies,  Ep  =  E\,  +  Ec.  For  general 
non-isobaric  configurations,  the  hot-spot  and  cold-fuel  pres¬ 
sures  are  linked  by  the  ratio:  pi,  =  <$>p  so  the  energy  of  the 
assembled  configuration  is  only  a  function  of  the  pressure: 

Etot  [ kJ ]  =  9.54  x  10 5pR3T3<t>-2pGl,r 

+  3.na°-6Mfmgp°Glr  (A3) 

In  conventional  ignition,  the  assembled  hot  spot  is  approxi¬ 
mately  isobaric,  since  compression  stops  and  the  pellet  begins 
to  expand  when  the  pressure  roughly  equilibrates  (<t>  ~  1).  In 
the  shock  ignition  case,  the  assembled  pellet  occurs  when  the 
pressure  in  the  hot  spot  is  higher  than  in  the  cold  fuel,  so  we 
expect  <J>  >  1 . 


FIG.  14.  For  the  isobaric  pellet  assembly  (<f>  =  1),  the  total  and  com¬ 
ponent  energies  of  the  pellet  assembly  as  a  function  of  assembly  pres¬ 
sure,  The  fuel  assembly  gain  (gray  curve,  with  axis  at  right)  is  also 
shown.  Note  that  the  peak  gain  occurs  at  a  pressure  that  is  greater 
than  that  at  the  point  of  minimum  total  energy. 


1.  Isobaric  or  shock  ignition  fuel  assembly 


At  low  pressures  the  hot  spot  energy  is  very  large,  and  the 
compression  energy  is  negligible  (Fig.  14).  As  the  pressure  in¬ 
creases,  the  hot  spot  energy  rapidly  decreases,  while  the  com¬ 
pression  energy  increases.  Conversely,  at  high  pressures  the 
compression  energy  is  dominant  and  increasing  with  pressure. 
This  implies  an  optimum  (minimum)  assembled  pellet  energy 
as  a  function  of  pressure,  which  is  found  at  the  point 


dE\ 

dp  )  Mf, a 


(A4) 


and  occurs  where 


Ec 

Eh 


=  5. 

isobaric 


(A5) 


The  pressure  and  energy  of  the  pellet  at  this  point  are: 

pmi"  [Gbar]  =  347  pRT5/4a~l/^~5/6Mj^g12  (A6) 

[kJ]  =  48  pRTl/2a1/2^l/3M5f^g  (A7) 

This  fully  determines  all  other  parameters  of  the  hot  spot  and 
cold  fuel  assembly  at  the  point  of  minimum  energy. 

The  point  of  minimum  pellet  energy  would  appear  to  give 
the  optimum  assembled  pellet.  However,  we  take  the  optimum 
point  to  be  at  the  maximum  pellet  gain,  which  is  near  (but  not 
at)  this  point.  (The  gain  here  is  defined  as  the  ratio  of  the 
yield  to  the  energy  of  the  assembled  configuration).  Assum¬ 
ing  a  burn  up  fraction  given  by  =  7/(7  +  pRtot ),  Fig.  14 

shows  that  the  point  of  maximum  gain  is  shifted  from  mini¬ 
mal  assembled  energy  to  a  higher  energy.  The  increased  yield 
due  to  a  larger  pR  more  than  compensates  for  the  increased 
energy  used  to  create  the  configuration. 


The  parameters  at  the  point  of  maximum  gain  can  be  found 
using  an  analysis  that  mirrors  Rosen’s27,  but  using  energy  in¬ 
stead  of  the  hot  spot  radius  as  the  independent  parameter.  We 
first  make  the  approximations:  (1)  £/,•  —  Ec  £),;  and  (2)  the 
burnup  fraction  can  be  approximated  by  the  simple  power  law, 
fbum  ~  '  (pR/7)1/2 .  Then  the  assembly  gain  is 

( oR)  i 

G  =  335.3  x  10 ifburnMmg/Eu  «  63366 Mmg^  J  (A8) 

Ekj 


With  the  first  assumption  we  can  relate  the  cold  fuel  pressure 
and  density  and  the  hot-spot  density  to  the  cold  fuel  energy: 


p  [Gbar]  =  0.036 a~3/2Mmg/2Ec/2 
pc  =  5 .45  (Ec/aMmg)3/2 

Pb^T^^2- 


Then  it  follows  that: 

Rh  =  r0E;5/1 
Rc  ~  r\Ecl/2  (Rl^Rl) 
pRrot  ~  5 .45  ( aMmg)~3/2Ec  x 
\n-r0E-2]+0.3pR 
G  ~  63366 Mmg  (^5A5(aMmg)~3^2E-1  x 

[n  -roE~2]  +0.3pR£/72)1/2 
where  we  have  defined  the  parameters 

r,  =  0.0 35ax/2MiJg  cm  —  ul/2 

and 

rQ  =  3].lpRTa3/2M5m/g2<t>'  cm  —  u5'2. 


(A9) 

(A10) 

(All) 


(A12) 


(A13) 

(A14) 
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We  have  also  assumed  that  the  cold  compressed  fuel  volume 
is  much  larger  than  the  hot-spot  in  the  equations  above. 

The  maximum  gain  occurs  at  dG/dE  =  0  while  holding  a, 
<f>,  and  Mmg  constant.  The  energy  at  maximum  gain  is  then 
E*  =  >/3ro/ri,  or: 

Etj  =  52  pRTl/2a}l2M5Jg<S>-112  (A15) 

where  we  have  assumed  pR  0.3 pR.  This  gives  the  maxi¬ 
mum  gain: 

G*  =  3150  pr1/4ilfi[4a'3/V/4  (A16) 

We  can  also  write  the  optimum  mass  as  a  function  of  energy 
by  solving  Eqn.  (A  15)  for  the  mass: 

M*ng  =  0.0087pW~3/5a~3/5<I>3/5E*f  ( A17 ) 

Other  parameters  of  the  pellet  at  the  peak  gain  point  are  found 
by  using  the  expressions  for  E^  or  M*  ;  they  are  listed  in 
tables  II  and  III.  We  note  one  curiosity:  at  the  point  of  maxi¬ 
mum  gain,  the  assembly  aspect  ratio  is  constant,  independent 
of  all  other  parameters:  Rc/Rh  =  '~^E*2  =  3. 

2.  Isochoric  fuel  assembly 

In  the  case  of  isochoric  assemblies,  p/t  =  p,  and  the  ratio  of 
pressures  between  the  hot  spot  and  cold  fuel  is  a  function  of 
the  cold  fuel  pressure  <t>  =  chip): 

«$>  =  153  Ta~3/5p~2/5. 

Equivalently,  we  can  write  this  ratio  as  a  function  of  the  cold 
energy: 

T>  =  575  TMmg/Ec. 

Because  of  this  functional  dependence  of  <f>  on  Ec,  the  min¬ 
imum  energy  changes  from  the  previous  result.  The  ratio  of 
the  energies  at  the  minimum  energy  point  becomes: 

Ec 

y  =3  (A18) 

isochoric 

which  matches  the  result  found  in  Eliezer  et  al.28.  The  pres¬ 
sure  and  energy  at  this  point  can  be  calculated  by  using 
Eqn.  (A3).  First,  rewrite  the  hot-spot  energy  equation  (Al) 
as  modified  by  the  relation  <J>  =  <!>(/?): 

Eh=  41  pfl3Ta6/V6/5.  (A  19) 

Then  the  cold  pressure  at  the  minimum  energy  is: 

pmm  =  8.8  pRl5/ST5/SMmg/8a3/ 8 
Using  Eqn.  (A2)  and  Etotai  =  ^Ec,  the  minimum  energy  is 

E^cJoric-mi"  [kJ]  =  12  pR3/4T1/4cc3/4M3/g4  (A20) 

Turning  our  attention  now  to  the  parameters  needed  for  op¬ 
timum  gain,  we  notice  that  the  expression  for  pRtoi  can  be 


simplified  in  the  isochoric  case:  pR,ot  =  p,  Rc-  The  radius  is 
still  given  by  eqn.  (A10),  so: 

pR,at  -  5-45  (aMmg)~3l2r\Ec  (A21) 

G  ~  27770  cT '/ 2M%3Ec/2/(Ec  +  Eh)  (A22) 

In  the  isochoric  assembly,  £),  can  no  longer  be  neglected 
compared  to  Ec  in  the  equation  for  gain.  The  gain  can  be 
simplified  if  we  express  £),  as  a  function  of  of  Ec  (using 
Eqn  (A  19)): 

Eh  =  2190  pR3Ta3M3ngE-3 

so  that 

G  ~  27770  a-l/2M2!gElc,2/(Ec  +  e 3E~3)  (A23) 

- 3 -  O  Q 

where  e 3  =  2190  pR  T  a  Mf  .  Then  the  maximum  gain  is 
found  from: 

dG  7  G  AGE3 

dE~c~  ~  2~Ec  ~  E*  +  e3 

which  has  the  root: 

E*c  =  (7e3)1/4 

or: 

E*  =  UpR3/4T1/4a3/4M3!g4. 

Note  that  at  this  point,  E%  =  ei/E*3  =  ^ E *,  so  that  EgJ  =  8£* 
at  maximum  gain,  or: 

E*kJ  =  12.7  pR3/4Tl/4a3/4M3lg4.  (A24) 

Invert  this  to  express  the  optimal  mass  in  terms  of  the  total 
energy: 

M*mg  =  0.0337  pR^T^a^E4^3.  (A25) 

Other  parameters  at  the  point  of  peak  gain  can  be  found  in 
tables  II  and  III. 


Appendix  B:  low  mode  uniformity  limits  and  convergence 
ratio 

We  consider  here  the  prospect  of  low  £  mode  perturbations 
distort  the  imploded  core  of  an  ICF  target  and  spoiling  gain. 
We  will  ignore  the  Rayleigh-Taylor  instability  (its  growth  rate 
is  ~  VI)  and  consider  only  the  secular  growth  of  initial  per¬ 
turbations.  Beginning  with  the  rocket  equation31  the  implo¬ 
sion  velocity  is  v  =  £ln[m(f)/mo],  where  m(t)  is  the  target 
areal  mass,  mo  is  the  initial  target  areal  mass,  P  is  the  driv¬ 
ing  pressure,  and  v  is  the  target  velocity.  The  implicit  as¬ 
sumption  is  that  P  and  the  mass  ablation  rate,  m,  are  con¬ 
stant,  so  m(t)  =  mo  —  mt,  or  t  =  =  (1  —  X) /X,  where 

X  =  m(t)/mo-  The  distance  pushed  is  then 

R~vt  =  —  \nX(\-X)/X2.  (Bl) 

mo 
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TABLE  II.  Fuel  assembly  parameters  at  peak  gain,  expressed  as  a  function  of  mass  Mmg 

parameter 

Isobaric  /  Shock  Ignition 

Isochoric 

c kJ 

52  pRTXI2all2M5Jg<t>-Xl2 

13  p3/4T1/4a3/4M3/4 

Ec/Eh 

21 /<b 

7 

p  [ Gbar\ 

705  'pRT5/4cc~l/4Mmg/U< L~5/4 

15p^15/8r5/8a3/8M~|/8 

Ph/P 

<L 

52  ^R~3/4T3/4a~3/4MlJg4 

Pc 

2037  pRT3/4a^^4M~l/4(i>~3/4 

202  p^9/8T3/8a-3/8Mffl3/8 

Phi  Pc 

0.09  pRV,1T~ 1  /2  a1  !2Mml  /6<J> 1  /2 

1 

Rh(pm) 

16  pRT  ]/4a'/4M5Jg'2<P1''4 

15p7T1/8T“3/8a3/8M3/8 

Rc(pm) 

3  Rh 

106  pTT^T^V^M11/24 

Gain 

3150  pRT~ 1  /4  a  "3/4M^4  3> 1  /4 

7280  p7?_3/8T_1/8a-7/8M3/g24 

TABLE  III.  Fuel  assembly  parameters  at  peak  gain,  in  terms  of  total  energy. 

parameter  Isobaric  /  Shock  Ignition  Isochoric 

Mmg  0.0087  pfif  0.034  pR~lT~1/3a~lE^3 

Ec/Eh  27/4>  7 

p[Gbar ]  5080  pRT3/2^~3/2E~l/2  132  pR5/2T5/6aEkJ5/6 

Ph/P  <$>  22p7T1r2/V1£1/ 3 

pc  6660  pf9/10ci"3/5r,/10£H3/10  120pR3/2Tl/2Eul/2 

ph/pc  0.2p"3/5r2/5a3/V/5£o1/5  1 

2.3  pf”1^1/2^2  4.2j^“1/2£^/2 

/?c(juw)  3fy,  22  p^~5/6r_Si/18a^1/3£'“/18 

Gain  963  p7^2/5q~9/10<f>2/5E3/10  2710  pR~2/3T~2/9 a-Vt’Ej/™ 


Different  drivers  give  different  mass  ablation  rates  and  mass 
remaining  at  the  time  of  stagnation.  Typical  values  for  final 
mass  are  Xfi,lai  ss  0.5  —  0.7  for  direct-drive  and  Xpnai  «  0.2 
for  indirect-drive  targets31 . 


1.  Sensitivity  to  initial  target  mass  asymmetry 


First  consider  asymmetry  due  to  surface  perturbations  or 
initial  areal  mass  perturbations.  (Thus,  inner-surface  or  outer- 
surface,  or  volume  perturbations,  are  considered  to  be  the 
same).  We  wish  to  find  the  variance  in  the  throw  distance, 
R,  as  a  function  of  mass  perturbation.  (This  variance,  when 
compared  to  the  assembled  target  size,  should  give  us  a  mea¬ 
sure  of  sensitivity).  Using  the  chain  rule  and  the  relations 
dX /dniQ  =  — X /mp  and  d X /dniQ  =  (1  —X)/mo,  we  find 


dR  _  R  1-X 
dniQ  mo  XlnX 


(B2) 


or: 


dR  X  —  1  dniQ 

R  XlnX  mo 


-MX) 


dniij 

mo 


(B3) 


X 


FIG.  15.  The  function  f\  =  ywc  describes  the  sensitivity  of  the 
core  radius  to  initial  mass  perturbations,  whereas  the  function  /2  = 
(1  +/l)/2  (in  gray)  describes  the  sensitivity  of  the  final  target  mass 
to  pressure  perturbations. 


The  function  fi(X)  relating  the  core  distortion  and  the  initial  The  allowable  initial  target  mass  perturbation  is  determined 
mass  asymmetry  is  shown  in  Fig.  15.  The  initial  perturbation  by  the  design  constraint  that  the  variation  in  push  distance 
dmo/mo  must  be  smaller  as  the  ablated  mass  increases.  (dR)  should  be  significantly  smaller  than  the  assembled  hot 
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spot  (Rh).  Defining  the  pellet  convergence  ratio  as  Cr  =  R/Rh, 
the  relevant  criteria  is: 


dR 

Rh 


CJ](X)—  <  1 
mo 


(B4) 


...for  some  meaning  of  “<C  1”.  Lindl32  notes  that  it  depends 
on  the  ignition  margin  available,  and  states  that  for  the  NIF 
ignition  capsules  it  means  “<  1  /4”.  This  is  roughly  consistent 
with  the  results  found  here.  So  we  are  left  with: 


dniQ 

mo 


1 

4Crfl(X) 


(B5) 


For  constant  density  shell  targets,  dmo/dm  =  dr/AR  where 
dr  is  the  surface  finish  and  A R  is  the  shell  thickness.  Typi¬ 
cal  convergence  ratios  for  ICF  targets  can  range  from  25-45, 
depending  upon  the  particular  target,  but  the  shock  ignition 
targets  here  have  larger  Cr. 


2.  Pressure  asymmetry  limits 


A  separate  issue  is  the  effect  of  low-mode  drive  pressure 
asymmetries.  We  can  repeat  the  previous  analysis  but  consider 
the  sensitivity  to  drive  pressure,  P.  We  use  eqn.  Bl,  and  also 
note  that  A  ~  y/P^ 1 .  Then  the  variance  of  the  distance  pushed 
with  respect  to  the  driving  pressure  is  given  by: 


dR 

~R 


1 

2 


i-x] 


X\nX 


dP 

~P 


(B6) 


The  function  /2(A)  is  shown  in  Fig. 15.  This  leads  to  the 
following  constraint  on  low-mode  pressure  variation: 


dP  ^  1 

P  <  4 Crf2(X) 


(B7) 


The  allowed  relative  pressure  variation  is  somewhat  greater 
than  the  mass  variation,  although  the  effect  is  only  appreciable 
as  X  =>  0. 
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